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Abstract 

We present calculations of nuclear level densities that are based upon the de- 
tailed microphysics of the interacting shell model yet are also computationally 
tractable. To do this, we combine in a novel fashion several previously dis- 
parate ideas from statistical spectroscopy, namely partitioning of the model 
space into subspaces, analytic calculations of moments up to fourth order di- 
rectly from the two-body interaction, and Zuker's binomial distribution. We 
get excellent agreement with full scale interacting shell model calculations. 
We also calculate "ab initio" the level densities for 29 Si and 57 Co and get 

reasonable agreement with experiment. 
PACS:21.10.-k, 21.10.Ma , 26.50,+x 
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Reliable nuclear level densities are important for the theoretical estimates of nuclear re- 
action rates in nucleosynthesis The neutron-capture cross sections are approximately 
proportional to the corresponding level densities around the nuclear resonance region. The 
competition between neutron-capture and (3 decay determines the fate of the 5 and r pro- 
cesses. Level densities can be extracted experimentally , but reaction network calculations 
require cross-sections for hundreds or thousands of nuclides, many of them short-lived. 

The most widely used description of the nuclear level density is the Bethe formula, based 
on a gas of free nucleon |§, and the modified "backshifted Bethe formula" f§]. Despite 
its ubiquity the Bethe formula is a phenomenological fit often requiring energy dependent 
parameters to match experimental data. 

The interacting shell model and other microscopic models accurately describe spectra and 
transition for a broad range of nuclides. On the other hand, "traditional" shell-model codes 
diagonalize the Hamiltonian in a large-dimensioned basis of occupation-state wavefunction 
but use the Lanczos algorithm to extract only a handful low-lying states. The level density 
requires complete diagonalization, a computationally forbidding requirement. 

An alternative to diagonalization is the Monte Carlo path- integral technique 0, which 
is well suited to thermal observables Although reasonably successful, path-integral 

methods are limited to interactions that are free of the 'sign problem' and are still very 
computationally intensive (i.e., require supercomputer time). Therefore we feel motivated 
to consider further alternatives. 

In this Letter we combine several previously disparate ideas based in nuclear statistical 
spectroscopy: (1) Analytic calculation of moments up to and including fourth-order; (2) 
partitioning the model space into subspaces; and (3) using binomial rather than Gaussian 
distributions. We find this combination to be successful, which we demonstate not only 
against exact shell-model calculations but also against experimental data. 

Nuclear statistical spectroscopy argues that many nuclear properties are controlled by 
low-lying moments of the Hamiltonian 0,0]. The first moment (centroid) is H — [L\ = \H\; 



2 



for n > 1 the nth central moment is p n = ^ (^H — H^j ^ , We also find it useful to introduce, 
for n > 2, scaled moments defined by m n = p n / {P2) n ^ 2 ■ The width p 2 provides a natural 
energy scale. 

Analytic formulae exist in the literature for computation of centroids through fourth 
moments directly from two-body matrix elements for any number of particles without the 
need for diagonalization in a many-body space [0-0]. Mon and French |fl showed that 
the level density in a finite space tends towards a Gaussian, which is described by the first 
and second moments. Any accurate description must however include deviations from a 
Gaussian, which require higher moments. The formulae for third and fourth moments fO 



are somewhat time-consuming (albeit less so than direct diagonalization and Monte Carlo 
path integration) and so, as far as we can divine from the literature, never implemented 
on a large scale. Grimes et al [JTJ] computed higher moments instead by the representative 
vector method, i.e., generating random sample wavefunctions in order to estimate the require 
averages. Today, on a modest workstation a few years old, we can compute, for a Ohu mid- 
pf-shell valence space, the third moments in a few hours and the fourth moments in a day 
or two (the time for a Ohu s<i-shell valence space is considerably faster). While not trivial, 
we emphasize such calculations, corresponding to roughly 10 8 levels, are still less demanding 
than Monte Carlo path integration by a factor of 10 or more in CPU time. 

Even with modern computers, however, moments beyond n = 4 are still numerically 
intractable. Another idea we borrow from statistical spectroscopy is partitioning the model 
space into subspaces, such as single-particle configurations (ld5/2) 4 (2si/2) 2 , etc. Then, rather 
than computing total moments as defined above, one calculates configuration moments. Let 
a denote a subspace and let P a be the projection operator for that subspace. Then the 
configuration moment is fi n (a) = (^P a (H — H^j ^ . Configuration moments are applied to 
configuration or partial densities, p a (E) = tvP a 8(E — H). The total level density p(E) = 
tr5(E — H) is simply the sum of the partial densities. 

Partitioning has several advantages. First, it comes at no cost: the formulae in the 
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literature for total moments are already expressed in terms of sums of configuration moments, 
which are thus needed and readily available. Second, higher total moments (n > 4) are 
dominated by lower configuration moments. Finally, partial densities themselves are useful 
for calculations of preequilibrium emission in compound nuclei [0.0 . 

The final step is a "base model" for the level density. Essentially, this base model 
makes assumptions for the values of the higher moments given the lower moments. The 
most common base model in statistical spectroscopy is the Gaussian distribution and its 
various modifications. (The random matrix model of Pluhar and Weidenmuller [15] is built 
upon semicircle distributions, which has both advantages and disadvantages that we do 
not have the space to discuss here.) A Gaussian is reasonable good starting point as it 
is already close to the actual level density. A common generalization is to expand in a 
Gram-Charlier series using Hermite polynomials 0; this unfortunately can lead to negative 
level densities. Another generalization is to extend the Gaussian to a function of the form 
exp(— aE 2 — j3E 3 — ^E 4 . . .) ||i~3| , |i~4"f but the relation between parameters a, f3, 7 and the 
moments is not amenable to a simple analytic formula. 

Instead we chose to follow Zuker, who recently gave a combinatorial argument that one 
should use binomial rather than Gaussian distributions to approximate level densities [|T7| . 
Zuker showed how third moments, that is, asymmetric distributions, are easily handled by 
binomials. Although Zuker did not comment on fourth moments, we find that the fourth 
moment can also be controlled in binomials. 

Consider the binomial expansion 

(i + a)" = £a».(^ (1) 

Now interpret this binomial expansion as the density of states. At the excitation energy 

E x = ek, e being an overall energy scale, the number of states is ^ k (^J ■ (Shortly we will 

see that deviations of A from 1 represent an asymmetric distribution.) Because we can write 

r^j with gamma functions, one can easily approximate it by a continuous distribution, 

n( E ) = \ Ex l e F(E max /e + 1) , . 

P{ x) r(E x /e + l)r((E max -E x )/e + l) {) 
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where E max = Ne represents exhaustion of the finite number of states. Although we began 
with N as an integer, it no longer has to be. 

The total number of states, which is the 'zeroth' moment, is d — (1 + A) , while the 
centroid and width are given respectively by 

NeX 

* = T+A (3) 
and the scaled third and fourth moments are 



m 3 = -7f=T=, ( 5 ) 



Vnx' 



For Gaussians (AT — > oo) m 4 = 3. For most binomials and for shell model diagonalization, 
the scaled fourth moment is less than 3, a typical value being around 2.8. In random matrix 
theory, the semicircle distribution typical of Gaussian Orthogonal Ensembles has m 4 = 2. 
Using Stirling's approximation, and a few others, Zuker arrives at 



/ in \ / ^ / /at i \ / -^a: i -^x . E max E x E max E x \ E x 

p ^ Ex) ~ v A^ exp r (iY_1) f — p — + ) +N F — lnA 

' J "l V J- 1 max ^max ^max ^max ' ^max 



(7) 

We carefully note that the above moments (0-0) are exact for discrete distributions but are 
only approximate for the continuous distribution (171). For large N, however, they are very 
good approximations. 

The key parameters of the binomial are the order N and the asymmetry parameter A. 
In the limit A = 1 and N —>■ oo one regains the Gaussian. If A = 1 then the binomial is 
symmetric: m 3 = 0; if A ^ 1 then the binomial is asymmetric. Zuker suggested that the 
order of the binomial, N, be fixed by the dimension of the model space. In that case N and 
A are fixed by solving d = (1 + A) and Eqn. (|5]) simultaneously. This we consider as the 
standard binomial, which can be asymmetric (nonzero third moment). 



We observe, however, that one could instead fix the order N by the fourth moment, 
and solve (||) and (||) simultaneously instead, afterwards multiplying the entire binomial 
distribution by a constant so as to get the correct total number of levels. This we refer to 
as the fourth-moment scaled (FMS) binomial. After N and A are determined, the centroid 
and width simply fix the absolute position and scale of the distribution. 

Elsewhere j0| we compare in detail the relative importance of third and fourth moments 
and configuration versus total moments. We find that the level density is best described by 
a sum of partitioned binomials, and that the difference between this best description and 
others is, at low energy, a factor of two or more. In the rest of this Letter we compare this 
statistical approach with exact shell model calculations, both from direct diagonalization and 
from Monte Carlo path integration, and with experimental data. Not only do we reproduce 
the secular behavior of the level density, we find that in some cases the partitioned model 
space can somewhat describe detailed structures at low excitation energy if one uses FMS 
binomials. 

To test the approach outlined above we considered a number of sd- and p /-shell nuclides, 
and show only a representative sample here; more can be found in |18|| . What we plot is 
in fact the state density, which includes 2 J + 1 degeneracies. Strictly speaking, the level 
density ignores 2 J + 1 degeneracies, but the literature is often cavalier with this distinction. 

First we consider full Ohu sd-shell calculations, using the Wildenthal USD interaction 
19| and focus on several nuclides that could be completely diagonalized using the OXBASH 



shell-model code IpOfl . The same single-particle energies and two-body interaction matrix 
elements were used by OXBASH for the exact calculation and the routines to compute the 
configuration moments for our statistical approximations. Fig. 1 shows three typical cases: 
32 S, 24 Mg and 22 Na. The histo exams are the exact state densities. We consider these three 
cases because they display very disparate collective behaviors: vibrational, rotational, and 
noncollective, respectively. While the standard binomial does reasonably well, the sum of 
partitioned FMS binomials describe all three somewhat better, despite different collective 
behaviors. We attribute the minimal differences for 22 Na to its noncollective, fully statistical 



behavior even at the lowest energies. The curious tracking by the summed FMS binomials 
of low-lying structure in the state density is very intriguing, although such tracking, as well 
as the structure in the exact histograms, is less pronounced in p/-shell nuclides. (NB: that 
all three examples are N = Z is irrelevant; the standard binomial and partitioned binomials 
work as well for N ^ Z although Gaussians fair even worse, because of more pronounce 
asymmetries (larger m^).) 

We have also performed Ohu calculations in the p/-shell. Figure 2 shows results for 
54 Fe and 48 Cr. Because of the prohibitively large dimensions of this model space, these 
nuclides cannot be diagonalized to yield all eigenvalues. Instead we turned to Monte Carlo 
path integration. To avoid the well-known sign problem || we fitted a schematic multipole- 
multipole interaction to the T = 1 matrix elements of the FPD6 interaction of Richter et. 
al |H] . Once again the sum of binomials is clearly superior to either a Gaussian or a single 



standard binomial. (It turns out that one needs a sum of binomials rather than a sum 
of Gaussians because the third configuration moments contribute nontrivially. Summing 
symmetric binomials also yields poor results.) We note that here the difference between the 
sum of "standard binomials" and sum of FMS binomials is negligible; furthermore there is 
much less low-lying structure in the state density than for sd-shell nuclides. 



Finally, in Figure 3 we compare directly to experimental data, for 29 Si |22| and 57 Co 
|23fl. Here it was necessary to include odd-parity states. For 29 Si we used the Wildenthal 



interaction [19] in the s<i-shell; to account for excitations out of the s<i-shell we included 



particle- hole excitations into the pf shell using the WBMB interaction of Ref. ||24|| . To 
avoid the u N-hu catastrophe" , we follow the suggestion of Ref. |24]| and decouple the nhuj 



excitations. Although clearly one should worry about contamination by spurious center-of- 
mass motion, we temporarily put this concern aside (as do all Monte Carlo calculations to 
date). Similarly, for 57 Co we used the modified KB3 interaction |25| and allowed 1-particle, 
1-hole excitations into a noninteracting g$/2 orbit whose single-particle energy is set by 
the start of abnormal parity states. Our calculations fall somewhat short at high energy, 
but this is clearly due to insufficient particle-hole states in our calculation. The overall 



agreement with experiment is remarkable, especially considering that the interactions were 
tuned to very low-lying states and not to more global properties such as state densities. 
While our treatment of particle-hole states here is admittedly ad hoc, keep in mind that 
nfiui excitations is always something of an art in shell model calculations. Indeed, we believe 
that our microscopic calculations are limited more by the uncertainty in the appropriate 
interaction than our statistical approximations, and that this should be the focus of future 
investigations. 

In summary: we have outlined a theoretical approach to level densities that is both 
microscopic in origin and also computationally tractable. Application to higher shells is 
hampered as much by our ignorance of the effective two-body interaction as anything else. 
In the near term we will also work to extend this approach to calculation of spin-cutoff 
factors and estimate of contamination by spurious states. 

This work was performed under the auspices of the Louisiana Board of Regents, contract 
number LEQSF(1999-02)-RD-A-06; and under the auspices of the U.S. Department of En- 
ergy through the University of California, Lawrence Livermore National Laboratory, under 
contract No. W-7405-Eng-48. 
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FIG. 1. Comparison of exact shell model state densities (histograms) from direct diagonal- 
ization in a full Ohu basis against statistical approximations for selected sd-shell nuclides. FMS = 
"fourth moment scaled" (see text). 
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FIG. 2. Comparison of exact shell model state densities (circles) from Monte Carlo path 
integration in a full Ohu; basis against statistical approximations for selected p/-shell nuclides. 
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FIG. 3. Comparison of experimental data of state densities for 57 Co and Si against a sum 
of binomials. Error bars in the lower figure represent the spread in experimental data. 
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